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Abstract. The dynamical justifications which lie at the basis of an effective Statisti- 
cal Mechanics for self gravitating systems are formulated, analyzing some among the 
well known obstacles thought to prevent a rigorous Statistical treatment. It is shown 
that N-body gravitational systems satisfy a strong chaos criterion, so supporting the 
assumption of an increasingly uniform spreading of orbits over the constant energy 
surface, i.e., the asymptotic evolution towards a microcanonical distribution. We 
then focus on the necessary conditions for the equivalence of statistical ensembles 
and remark that this equivalence is broken for any N-body system whose inter- 
action has a range comparable with its spatial extension. Once realized that the 
obstacles originate from the long range nature of the Newtonian interaction, leading 
to the non-extensivity of canonical and grand-canonical thermodynamic potentials, 
we show that instead a suitably generalized microcanonical ensemble constitutes an 
orthode, i.e., a reliable framework for a correct Thermodynamics. Within this setting 
we use the thermodynamic relations and find consistent definitions for Entropy 
(which turns out to be extensive, despite the non stable nature of the interaction), 
Temperature and Heat Capacity. Then, a Second Law-like criterion is used to select 
the hierarchy of secular equilibria describing, for any finite time, the macroscopic 
behaviour of self gravitating systems. 

Keywords: Gravitational Physics, N-body Systems, Statistical Mechanics 



1. Introduction 

Although is the most ancient among the known fundamental inter- 
actions, Gravity seems nevertheless to be the most reluctant to be 
included coherently within the theoretic frameworks describing the 
behaviour of microscopic and macroscopic many body systems, respec- 
tively. 

Indeed, in addition to the well known difficulties encountered by the 
attempts for a coherent formulation of a Quantum Theory of Gravity, 
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there are the less popular, though equally relevant from a concep- 
tual viewpoint, obstacles preventing a rigorous Statistical Mechanics 
of collisionless self gravitating N-body systems (1). 

This is due to the fact that the Hamiltonian function describing the 
gravitational interaction between N point particles of mass rrtj 

H(q, P ) = Y.^t ~ oEE I ' where r v = H " q<ll (!) 

i=i znii z i=i ^ n 3 

does not satisfy the stability property, which is a prerequisite for a self 
consistent Statistical Mechanical description and for the proofs of most 
of the theorems assuring the equivalence of the Statistical ensembles in 
the so-called Thermodynamic limit. The lack of ensemble equivalence 
for N-body self gravitating systems (SGS), leads naturally to the issue 
of determining which, if any, among the many ensembles complies with 
the orthodic condition (2). 

Notwithstanding the hindrances for a Statistical Mechanics of gravi- 
tational systems are known since many decades, many efforts have been 
devoted to overcome them, in the attempt to formulate, at least, an 
effective Statistical description of collisionless N-body SGS (3; 4; 5; 6). 

The astrophysical motivations for doing so are well evident looking 
at the otherwise inexplicable morphological, structural and dynamical 
similarities among hot stellar systems observed in the sky. The qualita- 
tive homogeneities of the luminosity curves and of velocity dispersion 
profiles, though with obvious quantitative differences, suggest unavoid- 
ably the existence of one or more relaxation processes leading globular 
clusters and elliptical galaxies towards an equilibrium state which is 
qualitatively independent on the global parameters of the systems and 
on the class of initial conditions of these astronomical objects at the 
formation era. 

Due to the peculiar nature of Gravity, the relaxation mechanisms 
(e.g., binary collisions) at work in usual laboratory many body sys- 
tems turn out to be unable to justify such a state of (quasi) equilibrium, 
because the corresponding relaxation time predicted is comparable to 
he Hubble time for globular clusters, and several orders of magnitude 
larger for elliptical galaxies. 

More than thirty years ago, Lynden-Bell, Saslaw and others (7; 8) 
first realized that the infinite range of Newtonian interaction could 
support relaxation mechanisms related to collective effects, in turn due 
to global, almost coherent, oscillations of the gravitational field, able 
to change appreciably the energies of the particles (i.e., the stars). 
However, the general relevance of such a mechanism has been revised 
and even criticized to some extent, (9; 10) and it is by now agreed that, 
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though useful to understand some aspects of the initial global evolution 
of gravitational systems, it alone cannot by sufficient to explain the 
universal light and velocity profiles observed. 

The lack of a satisfactory Statistical Mechanics for SGS made prob- 
lematic obviously any attempt to formulate a Thermodynamics for 
them. Among the various peculiarities of gravitational interaction (dis- 
cussed in the sections below), those implying the non extensivity of 
(dynamic and) thermodynamic potentials constitute the main obstacles 
to the definition of thermodynamic quantities and, even more, to derive 
standard relationships among them. 

A major piece of evidence of the very peculiar Thermodynamics 
proper to SGS is known since many decades to astronomers, but it 
went to know a wide popularity in the Physics community only when 
Thirring rediscovered this apparent paradox (11), which was already 
investigated rather deeply some years before by Antonov (12) and 
by Lynden-Bell and Wood (13). The problematic issue emerges when 
we recall from one side the well known (to astronomers) property of 
(bound) SGS of having negative heat capacities, and from the other 
the expression of the heat capacity as it can be easily derived within 
the Canonical ensemble (see, e.g., (14), §8.2): 

W + {{H -U) 2 )=0, (2) 

where U = (H) is the average value of the energy of the system in the 
canonical ensemble, f3 = 1/ksT is the (inverse) temperature parameter 
of the canonical distribution and ks is the Boltzmann constant. Written 
in terms of the heat capacity Cy, equation (2) reads 

C V = -Lj ((H 2 ) - (H) 2 ) > . (3) 

In the recent years, however, it has been realized that such a para- 
dox is not a peculiarity exclusive to SGS, but it is instead a common 
features of small systems, i.e., systems whose extension is smaller than 
(or comparable to) the range of the governing interaction (15; 16; 17). 
For these systems it is immediately realized that extensivity property 
breaks down. However these systems possess a well defined, though in 
some respects peculiar, thermodynamical behaviour, and so it turns 
out as well that such a property is not necessarily relevant to allow for 
a description of the macroscopic, quasi-thermodynamic behaviour of 
these systems, as discussed in more detail below, where we will analyze 
also why this feature obstructed to some extent the formulation of a 
Thermodynamic framework for SGS. 
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An immediate consequence of the lack of extensivity is the break- 
down of the ensemble equivalence, which is another milestone of the rig- 
orous statistical mechanical foundation of the thermodynamic formal- 
ism. However, the phenomenology of non extensive systems is currently 
investigated from theoretical, numerical and also experimental view- 
points, and many studies nowadays focuse on the dynamical, statistical 
and thermodynamic properties of Hamiltonian systems characterized 
by ensemble inequivalence (15; 17; 18). 

According to the considerations above, the attempt to formulate 
a coherent thermodynamic description of SGS appear to be feasible, 
starting from the basic assumptions which lie at the grounds of the 
Statistical description of many degrees of freedom dynamical systems, 
though paying due attention to every possible source of internal incon- 
sistency. 

Then, going back to Dynamics and Statistical Mechanics, we first 
investigate the requisites to be fulfilled by the Dynamics, in order to 
justify a statistical description of macroscopic properties of N-body sys- 
tems. This goal is accomplished using some concepts of Ergodic theory 
of (many degrees of freedom) dynamical systems, and in particular we 
adopt the Geometric Description of Dynamics, which allows to relate 
the chaotic properties of motions to the geometrical and topological 
properties of the dynamical manifolds associated to the system, see, 
e.g., (19; 20; 21; 22; 23; 24). 

Once found convincing indications that a self consistent effective Sta- 
tistical Mechanics of SGS can be formulated 1 , we then face the issue of 
the ensemble inequivalence and argue that the Classical Microcanonical 
ensemble alone is suitable for a reliable and coherent description of the 
Thermodynamics of non extensive systems like gravitational N-body 
ones. The criteria according to which the Microcanonical distribution 
over the hypersurface Eg of the full 6./V-dimensional phase space V, 
acquires such a privileged role over other classical ensembles are ana- 
lyzed and discussed, and it is shown how such a choice allows to derive 
a number of physically reasonable results on the macroscopic properties 
of N-body gravitating systems as well as unambiguous (re-)definitions 
of thermodynamic quantities and relations. 

As a by-product of the above analysis, we also derive and discuss 
some among the many non trivial relationships existing between the 
statistical properties 2 of dynamics over the constant energy hyper- 

1 Though with some peculiarities which are the statistical counterparts of the 
dynamical unique features of gravitational interaction. 

2 As, for example, ergodicity of dynamics in T space or even more strong features, 
like mixing. 
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surface Eg, the consequent statistical mechanical description, and the 
customary reduced description in terms of a distribution function (df) 
over the single particle six dimensional ^i-phase space, depending only, 
according to Jeans theorem (6), on the isolating integrals of motion 
admitted by the smoothed potential in the three dimensional configu- 
ration space. This theorem indeed is at the grounds of most of the works 
devoted to the construction of self-consistent models of stellar systems 
over the last four decades, and only recently it has been realized that 
its relevance for stellar dynamical modelling, at least in the customary 
strong formulation sketched above, has been somewhat overemphasized 
(25). We outline here the conditions to be satisfied by the full dynamics 
in FeAT-space in order to allow for the reduction process leading to a 
hydrodynamic (or kinetic) description in terms of a distribution func- 
tion over the /U-space. This represents a first step along the path of a 
critical reconsideration of the use of the single particle df as a tool to 
derive the equilibrium properties of stellar systems through maximum 
entropy like arguments based on it (26). A deeper discussion of this 
topic is however beyond the scope of this paper and will be the subject 
of forthcoming works. 

The paper is organised as follows: in section 2 we study the chaotic 
properties of many degrees of freedom systems, with particular em- 
phasis on the features of the dynamics of SGS and recall the usuals 
indicators adopted to quantify the degree of stochasticity of orbits in 
phase space. In order to relate Dynamics to the structure of the ambient 
space, we adopt the terminology and the tools of the Geometric Descrip- 
tion of Dynamics, which provides an elegant and effective framework 
to characterize the intrinsic properties of dynamical systems. 

To support our argumentations we recall some of the results ob- 
tained within this framework and reinterpret them in the perspective 
of the present work, devoted to propose thermodynamical description 
of SGS, justified on the grounds of the statistical features of dynamics 
and of the geometric properties of the dynamical manifolds. In sec- 
tion 3 we briefly introduce some concepts from the Ergodic theory of 
Dynamical Systems and use them to classify the statistical properties 
of the flow associated to the N-body dynamics, arguing for the pos- 
sibility of a statistical description of their macroscopic properties on 
any finite scale of times. In section 4 we review the issues related to 
the non equivalence of ensembles, then we discuss why the classical 
Micro canonical ensemble is the only legitimate to describe coherently 
the "Thermodynamics" of SGS; we then present some consequences 
of the thermodynamic relations derived, analyzing also the various 
inconsistencies encountered attempting to describe the system within 
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the canonical ensemble. Section 5 presents a brief discussion of the 
results, the peculiarities of the Thermodynamics of systems without 
Thermodynamic limit, a heuristic discussion on the concept of finite 
time equilibrium and concludes outlining the work in progress along 
the line here reported, possible extensions of the formalism and the 
improvements to the results presented which can be obtained with 
further studies and more extensive numerical work. 



2. Chaotic behaviour and Ergodic properties of Dynamics 
of gravitational N-body systems. 

Since the first rediscoveries 3 of the ubiquity of Chaos in nonlinear dy- 
namical systems in the 60's, (28), it was realized that unpredictability 
implied by the exponential sensitivity to small changes in initial condi- 
tions could be the ingredient needed to justify a more or less physical 
reformulation of the ergodic hypotesis, able to explain the approach to 
a stationary equilibrium state of most non integrable many degrees of 
freedom (MDoF) Dynamical Systems. To make a long story short, it 
is by now generally agreed that the presence of Chaos entails a system 
to forget the detailed memory of initial conditions, though some im- 
printing related to few global constraints or quasi conserved quantities 
could survive during the dynamical evolution, leading to a slow decay in 
time of specific correlation functions. The sufficient conditions implying 
irreversible behaviour and then able to justify a statistical mechan- 
ical description of macroscopic properties have been mathematically 
demonstrated only for very special (and rather abstract) classes of 
models, however there is a widespread agreement that, from a physical 
viewpoint, most MDoF nonlinear dynamical systems are indeed chaotic 
enough to relax at least to a stationary (meta-) equilibrium state. 

2.1. Chaos, Ergodicity and Mixing. 

For few dimensional systems Chaos is not sufficient to guarantee the 
efficiency of the approach to equilibrium, because diffusion in phase 
space can be inhibited by the fractal structure of deformed (or partially 
destroyed) invariant surfaces describing the integrable limit, i.e., by 
the Cantori originated by the disappearance of KAM-tori (29). Even 
above the KAM threshold, the exponentially long times associated to 

3 Most of the conceptual implications of dynamical instability as well as the 
realization that such a phenomenology should be the rule rather than the exception 
in the behavior of nonlinear (i.e., generic) dynamical systems were indeed already 
present in the Poincare work (27) more than a century ago! 
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the Nekhoroshev regime 4 (30) , the relaxation processes are practically 
frozen and from a physical viewpoint the approach to equilibrium will 
not take place at any finite time (31). 

However, both KAM and Nekhoroshev thresholds depends strongly 
on the dimensionality of phase space 5 , decreasing roughly exponentially 
with the number of degrees of freedom. Together with the Poincare- 
Fermi (32) theorem on the non existence of analytic integrals for non- 
linear systems with three or more degrees of freedom, the vanishing of 
the integrability or quasi integrability thresholds, suggests that MDoF 
chaotic systems should attain the asymptotic equilibrium states. 

However, from a physical viewpoint, it is also important that the 
approach to such states be fast enough, that is, it is required that 
the time averages of relevant dynamical macroscopic quantities con- 
verge to corresponding phase averages in finite times and not only 
asymptotically. Rephrased in terms of Statistical theory of Dynamical 
Systems, this amounts to say that dynamics leads any finite volume in 
phase space to spread out uniformly over the entire region accessible 
to motion. Formally, in addition to the ergodic property: 

limJ| T F[q(i),p(tp = F(q,p)//(dx) , (4) 

the dynamical evolution must also satisfy the stronger condition known 

as mixing: 

lim| r F[q(t),p(t)]G[q ) p] M ( ( lx) = Qf F(q, p)/i(*o) • (jf G(q, p)/i(dx) 

(5) 

where F and G are arbitrary functions of dynamical variables and fi(dx) 
is the invariant measure over T 6 . 

As above, the mixing property has been rigorously proved only 
for restricted classes of abstract dynamical systems, e.g., (29), though 
many theoretical and numerical works (33) now support strongly the 
belief that even more realistic models of many body physical systems 

4 During which the actions are almost exactly conserved and dynamical trajec- 
tories are trapped between confining surfaces. 

5 Throughout the paper we indicate with Af = N ■ d, the number of degrees 
of freedom of a dynamical system, N being the number of particles and d the 
dimensionality of the single particle configuration space. For Hamiltonian systems 
then the phase space T has dimension 2JV" and the constant energy hypersurface Eg 
has clearly dimension 27V — 1 . 

6 What is relevant from the point of view of Physics is that equation 5 be satisfied 
at least by macroscopic relevant functions F and not necessarily for any pair F and 
G, as instead is required in the mathematical definition of mixing. In the case of 
Hamiltonian systems the invariant measure is the Liouville measure. 
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possess good statistical properties and approach the equilibrium state 
on physically reasonable time scales; moreover, these statistical relax- 
ation times bear close correlations with the orbit diffusion rates in 
phase space, and also with Lyapunov exponents describing the degree 
of chaoticity. 

The above evidences have been given a formal expression recently 
through the formulation of the so-called Chaotic Hypotesis (34), which 
puts on the same level, for what concerns the statistical description 
of macroscopic properties, chaotic many degrees of freedom dynamical 
systems and Anosov systems, which are, among the abstract models, 
those characterized by the strongest statistical properties (together 
with Axiom A and K-systems, see, e.g., (35)) 7 . 

Before to proceed further, it is necessary to remark that all the 
above discussion about the KAM and Nekhoroshev thresholds is of 
fundamental relevance to single out the sources of Chaos in MDoF 
dynamical systems and for the understanding of the relationships be- 
tween the instability of the trajectories and the modifications occurring 
in the structure of phase space. Nevertheless, most of the tools and 
concepts developed within that framework are of little use for the 
characterization of the dynamical properties of SGS in an arbitrary 
state, because there is nothing indicating the existence of any integrable 
limit for gravitational N-body systems; i.e., for them, there is no hope 
to undertake a sensible perturbative approach to study the onset of 
Chaos, at least as long as In A" > 1. This amounts to say that the 
investigation of dynamical instability in SGS must resort to methods 
and tools pertaining more to the theory of strongly chaotic dynamical 
systems rather than to Perturbation Theory, whose approaches are 
suited for nearly integrable systems. 

2.2. Geometry of Dynamics and Chaos. 

Within the class of systems for which the mixing property has been 
rigorously proved there are the geodesic flows on manifolds of constant 
negative curvature (36), which were studied in the 30's by Hedlund 
and Hopf. The first to realize the relevance of such a line of research 
to Statistical Physics was Krylov (37), which studied the consequences 
of the instability of the geodesic flow associated to the dynamics of 
an hard sphere gas, predicting that the instability leads the system to 
macroscopically approach the equilibrium state, with a relaxation time 
which agrees rather well with more recent estimates. 

7 These models, though maximally chaotic, are, for what concerns the evolution 
of their statistical properties, almost as simple as integrable systems ((20), chapters 
4-=-6). 
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The formal apparatus for the geometric transcription of the dynam- 
ics of natural systems was developed independently by many authors 
and found a beatiful presentation in the Synge's article "On the Geom- 
etry of Dynamics" (38). Although completely unaware of the relevance 
of exponential instability, he there discussed extensively the stability 
of geodesic flows of dynamical interest. 

The fundamental argument able to relate geodesic instability and 
the entailed chaotic dynamics from one side, to the fulfillment of the 
mixing property on the other, goes through the proof that the dynam- 
ical manifold is characterised by almost everywhere negative sectional 
curvatures. If this is the case, then the Instability time, as measured by 
a suitable generalization of the (inverse) Lyapunov exponent, is closely 
related to the Relaxation time, giving the scale over which correlations 
decay and the memory of the initial state is lost. Then the point has 
been moved from a direct proof of mixing to the study of the sign of 
the sectional curvatures of the manifold. 

We refer for the details of the geometrical transcription of dynam- 
ics of Hamiltonian and Lagrangian systems to (24; 22) and references 
therein, and simply recall that the tool for investigating the stability 
of the geodesic flow are the Jacobi-Levi-Civita equations for geodesic 
deviation, which, in local coordinates over a suitable manifold read: 

Ts{^f) +n " AC = - ( a = 1, ... ,A0 (6) 

where <5q represents the perturbation to the flow, whose (possibly) 
exponential growth is the main signature of instability and depends on 
the stability tensor H. a c , defined along a given geodesic (with unitary 
tangent vector u), through the Riemann curvature tensor, R a bcd of the 
manifold (with metric g a b), by 

n a c = R a bcd u b u d ; (7) 

and V/ds is the total (covariant) derivative along the geodesic. From 
equation (6), an effective ordinary differential equation for the norm 
of the perturbation can be derived. Defining an unit vector v on the 
unitary tangent bundle of the manifold M, T q M, such that 5q = zu, it is 
straightforward, (22), to show that the magnitude of the perturbation, 
z, evolves according to: 

which is still an exact equation containing all the required informations 
about the overall instability of the flow. The simplification obtained 



Vv 



ds 



(8) 
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going from equation(6) to equation(8) is rather illusory. Indeed, even 
if we reduced a system of N differential equations containing covariant 
derivatives apparently to just one ordinary differential equation, the 
computational task has been made only a little easier. The problem is 
that equation (8) still contains the full stability tensor, whose explicit 
computation requires the knowledge of the 0{N 4l ) components of the 
Riemann tensor. In some particular cases it happens that the manifold 
is such that the simplification obtained is actually decisive (39) , but as 
expected, this is not true for most geometric transcriptions of realistic 
models of physical systems. 

For MDoF systems, however, in the hypotesis, to be verified a pos- 
teriori, that the dynamics is far enough from complete integrability 8 , 
it is possible to obtain a closure of equation(8), which then reduces to 
the approximate form: 

g+/C Ru [q( S )]z - 0; (9) 

where 

>c R M>)] = ^^; do) 

is nothing else than the Ricci curvature in the direction of the geodesic, 
equal also to 

R a bu a u b . Ric(u) 

where R ac = R b a bc = 9 bd Rabcd is the Ricci tensor. 

Referring to the above cited bibliography for a detailed discussion, 
we simply observe here how, from the effective equation (9) it is easy 
to understand the classical results on the chaotic properties of geodesic 
flows on manifolds of negative curvature: if the Ricci curvature per 
degree of freedom /CR u [q(s)] is everywhere negative along a geodesic, 
then the geodesic flow is exponentially unstable. So, a curvature every- 
where negative (that is, always, along a geodesic) entails an exponential 
instability of the flow. 



8 What is required actually is much less: it is sufficient that the evolution of the 
many degrees of freedom is not coherent, that is, that there is a spectrum of frequen- 
cies wide and far enough from resonances, so that collective periodic oscillations are 
avoided. In other terms, a good phase mixing is sufficient to prevent the persistence 
of phase correlations among the motions of individual particles; and we recall that 
phase mixing occurs in integrable systems with MDoF as well! From the geometrical 
viewpoint, this is enough to imply that the TV" — 1 non vanishing sectional curvatures 
are locally almost uncorrelated each other. 
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But this is not the only mechanism leading to Chaos. It is indeed 
possible to have an exponential growth of the solution of equation (9) 
also when the sign of the Ricci curvature along a geodesic is fluctuating 
and even for everywhere positive curvature. 

2.2.1. Sources of instability in positive curvature manifolds. 

Two further sources of instability are indeed possible for equations of 

the form (9): 

— the onset of exponentially unstable solutions can originate (when 
ICr > 0) from the mechanism of parametric resonance. This phe- 
nomenon, well known in the linear case (see, e.g., (40)), occurs 
for example in the generalization of Mathieu-like equations, 

x + n 2 (t) x = o, n 2 (t) = n 2 [i + f(t)} , </) t = o, (12) 

when the frequency S7 2 (i) undergoes a periodic or multiperiodic 
modulation whose frequencies satisfy suitable resonance conditions 
with the unperturbed average frequency, ft . If the behaviour of 
parametric resonance in linear systems is rather elementary, when 
this phenomenon couples with even small nonlinearities present in 
the system, the outcome can be very rich and complex (41). The 
interest in this mechanism resides in the fact that it appears very 
often when a phenomenological modellization of complex systems 
is attempted through a reduction process leading to the formula- 
tion of a system of effective equations, alike equation (9) above. 
It has been shown (e.g., (22; 23; 24)) that the occurrence of 
positive curvatures is the rule, rather than the exception, for most 
geometrical transcription of MDoF dynamical systems. As it is 
known, nearly all nonlinear systems exhibit nevertheless chaotic 
behaviours for almost all initial conditions. As every system is 
characterized by one or more intrinsic periodicities, often related 
to some kind of collective behaviour 9 , it appears reasonable to 
investigate the corrispondences between the parameter space and 
the astonishingly rich behaviour disclosed by evolutionary laws 
apparently so simple. In order to give a taste of the various phe- 
nomenologies displayed from an equation with only 1.5 degrees of 
freedom, we report in figure 1, the behaviour of different solutions 

9 As examples ranging from Solid State Physics to Astrophysics, we mention 
the resonances occurring with the lower modes frequencies in a slightly anharmonic 
crystal, which modulate the dynamics on the largest possible spatial scales, and the 
overall quasi periodicity of the gravitational field inside a stellar systems undergo- 
ing the Violent Relaxation phase, which modulate the background field of existing 
condensed substructures. 
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Figure 1. Parametric Resonance in a one degree of freedom nonlinear Lagrangian 
system with time-dependent potential. Each panel shows the plot of the solution of 
a single equation with a slightly different choice for the values of the only two 
free parameters. The behaviour of solutions of the first two panels recall those 
of the underlying linear model, with the typical damped and resonant evolution, 
respectively. Other choices of the parameters values originate either quasi periodic 
solutions, intermittent or steadily chaotic patterns. 

of a single equation, obtained adding small nonlinearities to (12) 
and whose two parameters have been varied slightly to obtain the 
different patterns shown in the panels. 

— The second mechanism able to lead to exponential instability of 
the flow is known as Stochastic Resonance and arises instead when 
the modulation of the frequency is due to stochastic fluctuations 
around its average value (42), and is typical of either non com- 
pletely isolated systems, where fluctuations arise from the envi- 
ronmental influence, or, as for parametric resonance, in the case 
of effective equations, describing the evolution of global or col- 
lective variables. Now the modulation is not described through 
deterministic modellization, but instead originates from the fact 
that the exact values of the parameters (and in particular of the 
average frequency) are determined by the combined effects of all 
the microscopic constituents of the system and are then influenced 
by the noise due to finite size effects. 

As a general rule (though of a rather heuristic nature), parametric 
resonance is responsible of the unstable behaviour of few dimensional 
systems as well as of the dynamical instability occurring in MDoF 
not very far from integrability and with an overall intrinsic underly- 
ing periodicity, whereas, as stated before, stochastic resonance is the 
typical mechanism promoting the development of Chaos in strongly 
non integrable MDoF systems; its description becomes more and more 
accurate as the number of degrees of freedom increases. 

It is obviously possible that the two mechanisms operate together, 
even if in this case it turns out very difficult to separate the effects 
and even to understand the conditions under which they enhance or 
damp each other. What is more plausible is that there are regimes in 
a given model in which either of the two mechanisms dominate over 
the other. If this is the case, heuristic argumentations suggest that in 
the intermediate range the two effects substantially do not interact, 
otherwise a discontinuous transition should be observed 10 . 

10 As the two mechanisms cause the onset of instability exciting in a very dif- 
ferent way the frequency space (strongly localized the first, with a rather narrow 
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From the above discussion we expect that the mechanism leading 
SGS to Chaos to be more closely related to Stochastic rather than 
parametric resonance, and we will see below that this is indeed the 
case. 



2.2.2. Geometric properties of the manifold and analytic computation 

of the maximal Lyapunov exponent. 
In recent years the geometrical description of dynamics has been ap- 
plied to the investigations of the dynamical behaviour of a large class of 
Hamiltonian and Lagrangian systems either with few or many degrees 
of freedom, (see, e.g. , (19; 20; 22; 23) and for a comprehensive review 
and an updated bibliography, (24)). 

The GeometroDynamical Approach (GDA) applied to Hamiltonian 
MDoF systems allowed to improve our understanding on many aspects 
of the characterization, quantification and sources of chaotic behaviour. 
Among the most remarkable, from the point of view of a concise syn- 
thetic characterization of the dynamic instability, it has been possible to 
use the GDA, in conjunction with the theory of Stochastic Differential 
Equations (42), to derive an effective analytic formula for the maximal 
Lyapunov exponent (21), which is obtained simply evaluating average 
values, fluctuations and autocorrelation time of the frequency of the 
effective equation (9). 

If we indicate with Ur = (/Cr u [q(s)]) the average (positive) value 
of the Ricci curvature per degree of freedom, with o\ r = M(t\ r the 
rescaled variance of its fluctuations 11 and with t^ r the autocorrelation 
time of the stochastic process describing the time evolution of JCr u [q(s)] 
along any given geodesic, it is possible to analytically compute the 
average exponential growth rate of the perturbation according to the 
following formula 



where 




^ = \F? U" ' . < 13 > 



7] 



and 



9 a l R T k R 



(15) 



distribution the second), we expect the effect of the overlap to be in generic cases 
negligible. 

11 For a discussion of the M rescaling of the variance, see (21). 
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The importance of equation (13) relies also on the fact that for most 
many body systems the quantities entering it can be computed using 
either the standard time averages along numerical trajectories or the 
statistical phase averages over the canonical or micro canonical distri- 
bution. The agreement between the results obtained using either of 
the above procedures is rather satisfactory (20; 21), and numerical 
simulations along trajectories allow also to check the reliability of the 
assumptions at the grounds of the use of the theory of stochastic dif- 
ferential equations, first of all the gaussian distribution of the values of 
the curvature. 

In addition to this, the GDA helps to improve the conceptual un- 
derstanding of the onset of instability in MDoF Hamiltonian systems, 
showing that, even for the most chaotic N-body systems, the sources of 
instability rely mostly on the stochastic fluctuations of the curvature 
rather than on the occurrence of negative values. In particular, for 
what concerns SGS, it has become clear that dynamical instability of 
trajectories, originates mainly from binary encounters, and that the 
latter entail positive fluctuations of the curvature: the strong chaoticity 
of N-body systems is mainly due just to the huge amplitude of cur- 
vature's fluctuations around the average value, to be compared with 
what happens in other MDoF Hamiltonian systems, where fluctuations 
become of the order of the average values only in the strongly chaotic 
limit. 

Moreover, the GDA singles out clearly the different regimes accom- 
panying the approach to equilibrium, as the transition to a less unstable 
regime when virial equilibrium has been attained and the violent re- 
laxation process accomplished; and many other phenomena operating 
on longer timescales and finer spatial extensions. The details of these 
issue have been partially discussed in (22) and will be analysed in a 
deeper detail elsewhere. 

2.3. Strong Stochasticity threshold and Chaotic 
Hypotesis. 

To the main issue of this paper what is relevant of the above discussion 
is the evidence of a very strong chaoticity in gravitational N-body 
systems. Without going into formal details, we present here the results 
of the application of the above formalism to SGS, and compare them 
with analogous outcomes related to a chain of anharmonic oscillators 12 , 
in order to characterize the degree of chaoticity present in N-body 
dynamics. As anticipated above, in MDoF nonlinear systems there are 
in general various dynamical regimes. Starting from quasi integrable 

12 The well known Fermi-Pasta-Ulam quartic, or FPU-/3, model (32). 
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behaviour for small nonlinearity (or at enough low energy), a complete 
hierarchy of chaotic regimes is encountered as the departure from in- 
tegr ability increases. The KAM and Nekhoroshev thresholds are very 
near the zero of the perturbation parameter, e, as soon as the number of 
degrees of freedom is greater than a few tens, but nevertheless, in sev- 
eral different Hamiltonian systems it is possible to single out a further 
physically meaningful threshold, separating two different regimes of 
stochasticity. This dynamical transition has its statistical counterpart 
in a rather abrupt change in the behaviour of the relaxation time of 
macroscopic quantities. 

Briefly, below such a threshold, let it be e c , the rate of increase of 
the maximal Lyapunov exponent, Ai(e), is described by a relatively 
steep (model dependent) power law and the relaxation times are very 
long, depending sensitively on the perturbation parameter. Above this 
threshold, the growth rate of Ai (e) slows down and the relaxation times 
become much shorter and nearly independent of the value of e (as long 
as £ > e c ) (33). 

A very interesting and illuminating feature of the scaling law Ai(e) 
above e c 13 is such that the dimensionless quantity 71 = Ai • to, intro- 
duced in (22), turns out to be virtually independent on the perturbation 
parameter. This because, above the SST, the maximal Lyapunov expo- 
nent scales exactly as the reciprocal of the dynamical time scale. This 
evidence implies that the diffusion of orbits in phase space proceeds at 
the fastest rate allowed by the dynamics. This fact has been noticed and 
briefly discussed in (22), suggests a very elegant scenario for the mech- 
anisms governing the onset and the development of Chaos, and gives 
a convincing picture of the modifications in the phase space structure 
accompanying the departure from integrability. The investigation of 
these issues and of their consequences is however beyond the scope of 
these pages, and will be discussed elsewhere. 

To proceed towards our goal, i.e., a Thermodynamics for SGS, in the 
figures 2 and 3, we compare the energy dependence 14 of the maximal 
Lyapunov exponent, computed either using the standard algorithm (43) 
or the analytic geometrodynamic formula, equation (13). 

Figure 2 displays the transition in the scaling behaviour of 71 (/3e) 
outlined above. 

In figure 3 are reported the values of the maximal Lyapunov expo- 
nent as a function of the specific energy of the system, for two series of 
simulations with N = 400 and N = 4000 bodies. It is evident the power- 

13 The so called strong stochasticity threshold, SST, (19; 33). 

14 To be precise, for the FPU model, the perturbation parameter is the product 
(3e, where /3 is the anharmonicity constant. 
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Figure 2. Scaling of the dimensionless chaoticity indicator 71 = \\to with the 
perturbation parameter fje for a FPU chain, as obtained according to the analytic 
formula, eq.(13). Open circles refer to the values computed from mean values and 
fluctuations obtained from phase space averages (21), whereas, starlike symbols with 
error bars show the values obtained from time averages computed along numerically 
integrated trajectories (20). 

Figure 3. Scaling with the binding energy per particle, e, of the standard maximal 
Lyapunov exponent Ai and of the dimensionless chaoticity indicator 71 = Ai to for 
two sets of gravitational N-body systems, with N = 400 (circles) and N = 4000 
(stars) bodies, respectively. 

law dependence of Ai(e) oc |e| 3 / 2 . This dependence however disappears 
if we measure the instability exponent in the unit of the dynamical 
time to, which is customarily defined for SGS as proportional to the 
inverse square root of the density try oc (Gp)^ 1 / 2 . So, for SGS, we find 
that the adimensional chaoticity parameter 71 is constant and does not 
depends on the amount of binding energy in the system. 

From the comparison of the two figures the peculiarity of SGS emerges 
clearly: they behave, at any binding (specific) energy, as standard many 
body systems above the strong stochasticity threshold. 

This not only confirms that bound collisionless SGS are character- 
ized by a dynamical instability proceeding at a very fast rate 15 , but 
also suggests that they possess the strongest statistical properties, in 
analogy with those of standard dynamical systems in the regime of fully 
developed stochasticity. 

Once realized this, we are allowed to invoke the Chaotic hypotesis 
outlined in the Introduction, and to build a Statistical Mechanics for 
SGS. 



3. Chaos and Statistical properties of Dynamics. 

If the dynamics of a MDoF system possess strong enough stochastic 
properties, the Chaotic hypotesis (34) asserts that the averages of phase 
space functions can be computed as if the system were a transitive 
Anosov flow. From a physical viewpoint, this amounts to say that 
the system attains the most probable macroscopic equilibrium state 
compatible with the global constraints imposed at the initial time 16 . 

15 This is indeed a very well known fact: it was discovered already by Miller in 
the '60s, when he performed the first pioneering numerical simulations of N-body 
gravitational systems and has been confirmed many times also in the last years. 

16 These constraints are essentially the values of the possible conserved quantities, 
like energy, linear and angular momentum and so on. 
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So, if we are able to find this equilibrium state, we can compute all 
the averages using the invariant measure over the stationary distribu- 
tion. The analogy with Anosov systems says even more, it gives also an 
indication about the rates governing the approach to the equilibrium 
state, which should be exponential if the manifold is compact. 

For most MDoF systems these equilibrium (or stationary non equi- 
librium) states coincide either with the canonical or microcanonical 
ensemble, according to the system being in contact with a thermal 
bath or isolated, respectively. If the interaction potential governing the 
internal dynamics satisfies the conditions outlined in the Introduction, 
then the theorems of Statistical Mechanics entail that the descriptions 
based on the two ensembles agree, apart from terms which vanish as 
1/N in the thermodynamic limit, i.e., when the number of particles, the 
spatial volume and the energy of the system grow indefinitely, keeping 
however constant their ratios: 

N ^ oo ; F ^ oo ; £ ^ oo (16) 

with 

■ N ■ ? 

n = — = constant ; e = — = constant . (17) 

3.1. Obstacles to a Statistical Mechanics for gravitating 

SYSTEMS. 

From the above discussion, one amongst the many difficulties for a Sta- 
tistical Mechanics of gravitationally bound SGS emerges immediately: 
it is not possible to perform the thermodynamic limit complying with 
the above requirements. 

From the expression of the gravitational N-body Hamiltonian, equa- 
tion (1), it is immediate to realize that keeping constant the density 
n implies that the specific energy grows without limit, e oc N 2 / 3 . This 
problem derives obviously from the fact that, for bound SGS, the num- 
ber of bodies (or better, the total mass), the energy, and the spatial 
volume cannot be fixed independently from each other, as it is instead 
possible for standard laboratory systems. 

It is usually stated that that the one just mentioned is only the most 
evident difficulty and that another one, in principle not less relevant, 
is related to the fact that the most probable state for the gravitational 
N-body system cannot be determined, because formally, the measure 
of the phase space volume accessible for a given total energy is virtually 
infinite. 

However, if we are not interested in the properties of the abstract, 
eternal, final equilibrium state, but rather to determine the quasi sta- 
tionary states whose secular evolution proceeds on time scales much 
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longer than any dynamical process , the second of the above mentioned 
difficulties can be surmounted, as we will show, in a (hopefully) self 
consistent way. 

The remarkable fact is that, to overcome this obstacle, we do not 
need to avoid or mask the peculiarities originating it, but rather we 
succed to incorporate their implications, to outline a hierarchical evo- 
lutive path, able to describe the succession of quasi-equilibria of stellar 
systems, covering a range of timescales spanning from those of the order 
of the virialization (or Violent Relaxation) up to those related to the 
formation of a collapsed-core-expanding-halo structure. As long as this 
hierarchy of timescales is covered, dynamical processes involving finer 
and finer spatial scales take place and this evolution is actually endless, 
as the final equilibrium state can be attained only on infinite times, and 
just for this, does not bear any direct relevance to real astronomical 
objects, once verified that the secular evolution of the metastable states 
proceeds on timescales longer than the stellar systems lifetime. 

From the above discussion, it follows that the sole serious diffi- 
culty for a standard statistical description of N-body SGS is related 
to the non existence of the standard Thermodynamic limit. As we 
will discuss below, this is not a real hindrance for a self consistent 
Thermodynamics of SGS, as it is well known that many other labo- 
ratory systems 18 , characterized by non-extensive interactions, do not 
possess a well defined Thermodynamic limit 19 and, nevertheless their 
macroscopic or collective behaviour can be described coherently using 
a suitable thermodynamic formalism which take into account their 
peculiarities. 

The analysis below shows that a conceptually coherent framework 
can be settled to describe also SGS, certainly intrinsically even more pe- 
culiar than the small systems outlined above, but nevertheless suitable 
for a self consistent Thermodynamics. 

3.2. Dynamical singularities and Statistical Mechanics. 

The above scenario owns also a dynamical counterpart which can be 
easily understood within the Geometro dynamical framework. With- 

17 And, from the astrophysical viewpoint, also much longer than the lifetimes 
of stellar systems to which we plan to apply at the end all these theoretical 
considerations. 

18 The so called "small systems", see below. 

19 Rather, the point is that, loosely speaking, the behaviour predicted in the T.L. 
does not bear any resemblance with the phenomenology observed in real, finite size, 
systems. And this should not to be considered surprising, because, for those systems, 
some peculiarities come just from the finiteness and disappear when the size exceeds 
the range of the interaction. 
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out going too much into details, we remark that the convergence of 
time averages to phase space values of geometric quantities (e.g., the 
curvature) proceeds in the case of SGS on timescales which strongly 
depend on the spatial resolution over which we compare the averages 
themselves. That is, the time required for the time averages (and fluctu- 
ations around them) to converge to the static averages is very sensitive 
to the spatial resolution used to compute the latter. Incidentally, we 
remark that the ergodicity time for the Ricci curvature is infinite if 
the mathematical definition of the Newtonian interaction is assumed. 
It becomes viceversa a well defined, finite, quantity, if the physical form 
of the interaction is considered. It must be stressed that even in this 
respect, as the spatial resolution increases, the ergodicity timescales of 
geometrodynamic quantities grow up, exactly at the same rate of the 
times associated to the relaxation of statistical quantities. 

3.2.1. Mathematical vs. Physical peculiarities of gravitational 
interaction. 

We conclude this section listing the more relevant peculiarites of the 
Newtonian interaction outlining also which, if any, amongst them could 
constitute a real obstacle to a functional statistical mechanical founda- 
tion of a Thermodynamics for SGS. 

Singular interaction at short distances: The most evident singularity of 
the gravitational interaction is obviously related to the divergency 
of the potential when two particle get closer and closer. Although 
clearly of outmost importance for the single particle dynamics and 
for the evolution of few body systems, the global and statistical 
relevance of close encounters is practically negligible over any finite 
timescale, as the frequency of high energetic close encounters is 
vanishingly small and the effect on the global quantities irrelevant, 
compared with that of the overall dynamics. Then, even if it is 
physically plausible 20 , from a statistical viewpoint there is no need 
of any softening of the potential. 

Non confining nature of the potential and divergence of the phase-space 
volume: The configuration space of a N-body bound SGS obviously 
diverges, and in principle very badly, indeed as 

r p ~ lim v N - 2 . 

This is is due to the fact that Newtonian interaction is unable 
to confine, asymptotically, inside a finite volume more than two 

20 And numerically very useful, especially if low order efficient symplectic inte- 
grators are used, in order to guarantee the faithful preservation of the Hamiltonian 
properties of the flow. 
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bodies. If the mathematical N-body system is considered, an anal- 
ogous divergence comes out also from the volume in the momentum 
space. However, as long as logiV S> 1, the characteristic time 
for the evaporation is considerably longer than any dynamical 
timescale, and the evolution is only secularly influenced by the 
slow evaporation of stars, which, incidentally, causes the remaining 
system being more and more bound, so inhibiting further escapes. 

Unscreened, long-range nature of the interaction: It is the singularity 
which most affects the statistical and thermodynamics proper- 
ties of SGS. This is in part obvious as it concerns the dynamics 
over all the spatial scales, and consequently, according to the gen- 
eral discussion above, every timescale. The long range nature of 
the interaction reflects on the efficiency of collective phenomena, 
among which the Violent Relaxation is surely the most popular 
but it is not the only one. The absence of any effective screening 
mechanism entails that gravitationally bound systems do not own 
any characteristic lenght, and reflects on the scale invariance of 
most features. From statistical and thermodynamical viewpoints 
the most far reaching consequences reflect in the following points: 

— N-body systems are not stable in the sense of Rigorous Sta- 
tistical Mechanics (1): given any number of particles N, the 
binding energy per particle can increase without any upper 
bound (i.e. the energy does not have a lower bound). In 
particular, the binding energy of the fundamental state of 
a system scales at least as N 2 . That the lack of stability 
property is related to the long range nature of the potential 
and it is unaffected from the singularity at short distances 
emerges clearly from the discussion above and from many 
numerical simulations performed recently (44). Indeed, the 
introduction of a softening of the potential is able to introduce 
an upper limit for the binding energy of the system, but is 
unable to cure the lack of the stability property 21 . 

— From the previous point, it follows that the internal energy of 
a system scales as N 2 , and this remains true for most thermo- 
dynamic potentials. As the extensivity of the latter is a neces- 
sary condition for the proof of the ensemble equivalence, the 
long range nature of the Newtonian interaction is the main 
cause of the peculiarity of Gravitational Thermodynamics. 

21 Even the most radical modification of the short range interaction, i.e., the 
introduction of an hard-core repulsive potential at short distances, does not prevent 
the divergence of the binding energy per particle in the thermodynamic limit, though 
being able to make it finite at any finite N. 



CPMPfinal.tex; 1/02/2008; 17:02; p. 20 



Prom Chaotic Dynamics to Thermodynamics of self gravitating systems. 21 



From the above discussion it emerges clearly that is the non-extensivity 
of potential energy (and, by the virial theorem, also of the kinetic and 
total energies) to imply the non equivalence of statistical ensembles, 
and this requires, in turn, a deeper analysis in order to understand 
which ensemble, if any, is best suited to describe SGS. 



4. From Statistical Mechanics to Microcanonical 
Thermodynamics of self gravitating systems. 

Up to now we worried about the rigorous theorems concerning the 
equivalence of ensembles, the scaling of macroscopic quantities in the 
thermodynamic limit and so on. However, as emphasized in previ- 
ous sections, in recent years it has realized that the issue of the in- 
equivalence of statistical ensembles emerges whenever the scale of the 
dominating interaction is comparable with the spatial extension of the 
system under study. It is easy to understand that in such circumstances 
the interaction energy ceases to be an extensive quantity and, if some 
equipartion property holds, the same is true also for kinetic and total 
energies. However, many of these so-called small systems (15) possess 
a well defined, though somewhat peculiar, thermodynamic behaviour. 
In these systems it is not difficult to encounter regions in the space of 
thermodynamic parameters where the specific heats become negative. 
When described within the microcanonical ensemble, all these peculiar 
phenomena do not originate any paradox, which appears instead as 
soon as a canonical ensemble approach is attempted. The best that 
can be done, sometimes, within this framework, is to describe the 
phenomenology inside these regions as a phase transition, though true 
phase transitions themselves are rigorously defined only in the Ther- 
modynamic limit, which we know to be a rather meaningless limit just 
when the system properties are not extensive. 

4.1. Ensemble inequivalence. 

As discussed in the previous sections, the non extensive nature of the 
energy entails, among others, the inequivalence of statistical ensembles. 
In order to find out what among the classical ensembles is best suited 
to describe SGS we first observe that, if we consider isolated N-body 
systems, nothing suggests to suppose them in contact with any heat 
bath. Going a little more below the surface, we observe more concretely 
that if we require, according to Boltzmann, that the ensemble fulfils 
the orthodic property (see (2), §5), then neither the canonical nor 
the grand-canonical ensembles can be used to describe gravitationally 
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bound N-body systems. Heuristically, this simply because both these 
ensembles fix the Temperature as a constant parameter of the distri- 
bution, assuming implicitly that this quantity is constant throughout 
every system of the ensemble. 

To this somewhat formal point, we have to add the already quoted, 
and physically most relevant, incongruence met by the canonical (and 
also grand-canonical) approach: SGS are characterized by negative heat 
capacities, whereas equation (3) shows that within the canonical en- 
semble the heat capacity is positive by definition. Moreover, the non 
extensivity of potential, kinetic and total energies entails the analogous 
non extensivity of canonical and grand-canonical thermodynamics po- 
tentials, as the Helmholtz free energy F or the Gibbs potential G. 
These features not only imply the non equivalence of ensembles but 
also indicate that even the definition of averages (and even more of 
fluctuations) of thermodynamic variables becomes problematic. 

Up to this point we analyzed the problems affecting the canonical 
and grandcanonical descriptions. Now we turn to ask how the pecu- 
liarities of gravitational interaction reflect also on the microcanonical 
description. 

Before to proceed, we think that this seems to be the appropriate 
point to recall a general statement: the strong chaotic properties of 
dynamics, the Poincare- Fermi theorem, the general theorems of Ergodic 
Theory, all focus on the relationships between the stochastic proper- 
ties of motions and the occurrence of a uniform diffusion of initially 
concentrated phase space volumes over the whole region accessible to 
motion, compatibly with the integrals of motion allowed by the system 
dynamics. If such a spreading process in phase space is efficient enough, 
then, those general theorems, also according to the Chaotic Hypotesis, 
assert that the the time averages of dynamical variables coincide with 
the phase averages computed over the region of phase space where 
motions take place transitively. If we are dealing with Hamiltonian sys- 
tems, this region is the constant energy surface 22 Eg. This amounts to 
say that the strong enough stochastic properties of Dynamics of MDoF 
Hamiltonians, naturally lead to a microcanonical statistical description. 

The task of the powerful formal apparatus of rigorous Statistical 
Mechanics is then to prove the equivalence of ensembles. In case of 
success, then we can use the more tractable expression of canonical 



Possibly taking also into account the conservation of other global integrals, like 
the position and linear momentum of the center of mass of the system and the total 
angular momentum if only central forces are present. This reduces the dimensionality 
of Sf from 6iV — 1 to 6iV — 10, which, for log N ^> 1, does not modify obviously the 
macroscopic behaviour. 
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or grand-canonical ensembles to describe the macroscopic properties 
and the Thermodynamics of our system. If the equivalence does not 
hold, then the sole Statistical description allowed is the one compatible 
with the global constraints on our system: if it is isolated, then only a 
microcanonical description is legitimate. 

From the above discussion, it is clear that the use of microcanon- 
ical description is justified, provided that no internal inconsistencies 
occur; and that the inequivalence of ensembles, if ever, should suggest 
a critical reconsideration of the results obtained using a canonical or 
grand-canonical description. 

It is easy to see that the only trouble comes from the formal un- 
boundness of the phase space volume accessible at a given constant 
energy. Nevertheless, if we are interested in the finite time quasi equilib- 
rium states, we can take the argumentations discussed in the previous 
sections to emphasize that, at any finite time the effective phase space 
volume accessible is well defined in terms of the global quantities de- 
scribing the state of the system (energy per particle, specific angular 
momentum, etc.) and it is finite 23 . 

So, for any finite scale of time, we are able to derive an effective 
expression for the microcanonical entropy, secularly evolving, and able 
to take into account the increasing effects of the phenomena occurring 
on smaller and smaller spatial scales, and, correspondingly, on longer 
timescales. 

Starting from the Hamiltonian of the system (1), given a value for 
the total energy of a bound SGS 24 , £ < 0, we can define a set of 
macroscopically averaged quantities: 



the energy per particle: 



e = | < , (18) 



It is easy to verify that even in the presence of escaping particles or taking 
into account the occurrence of close encounters, the contribution of these events 
to the actual phase-space volume is negligible as long as log N 3> 1 for times of 
t>0(Nt D ). 

24 It is possible to perform a statistical mechanical treatment even for loosely 
bound systems or systems in interaction with some background gravitational field. 
This leads to a very interesting phenomenology, indicating the occurrence of critical 
phenomena in 3D-gravitational systems. The details of this work are described in 
(44). 
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the virial ratio: 



Q(t) = 



-2K{t) 



N 2 

v £■ 

i=l 



N 



(19) 



1 V m ' m J 

2 2^ 2^ r . 

i=\j+i 13 



where the indication of the time dependence of the kinetic and 
potential energies, K(t) and U(t), respectively, is a shortcut to 
recall that their values depend on the state of the system at time 
t; 



the inertial radius: 

Ri{t) = 

the harmonic (or virial) radius: 

R g (t) = 



N 

E m i r i 
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N 

E m i 
i=l 

N 

EE™ 

i=ijjti 



1/2 



(20) 



(21) 



and the ratio among them: 



a(t) 



Ri(t) 

Rg(t) 



(22) 



In the simplifying hypotesis that the mass spectrum is nearly flat, 
so that the average mass in = (mj)j represents the typical mass of the 
system 25 , it is a simple computation to obtain an expression for the 
microcanonical entropy of the system, which reads 



S = ks ct(N, e, Q, a, m) , 



where 



3 AT 
a=— In 



G 2 N 2 m 5 
-2s 



a\t)Q{t) [2-Q(t)\ 



(23) 



(24) 



5 When a very steep spectrum of masses is allowed, the computations are only a 
little more involved, but the results presented below are unchanged. In addition it 
is possible to obtain also the dependence of the entropy on the mass sectrum. 
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and it is obviously understood that the also the varying quantities a(t) 
and Q(t) depends on time only through the actual state of the system. 

The above argument can be extended to take into account other de- 
tails about the state of the system, but the simplest derivation sketched 
here is sufficient to emphasize a number of relevant points. 

It is evident that the microcanonical entropy (24) gives correctly the 
equilibrium conditions related to variations of the parameters on which 
it depends. 

However, the most striking result is that it turns out to be an 
extensive quantity, being proportional to the number of particles 26 . 

Another very important feature of this entropy is that it is possi- 
ble to derive from it easily all the relevant thermodynamic quantities, 
which agree with the phenomenological expressions obtained at a rather 
heuristic level. 

The most important is the microcanonical temperature which, co- 
herently with the statistical description, is clearly only an average 
quantity for the systems of the microcanonical ensemble 27 

1 . fdS\ 1 fdS\ 



T \d£j N N\de /N 

which, using (24), gives, coherently with the fact that bound systems 
have negative energy and positive temperature, 

e = -h B T . (26) 
Consequently, a negative heat capacity follows: 

^ (If) „ = (27> 

that is, a negative specific heat, cy = —3/2ks- 

The surprising features of the entropy (24), defined in such a straight- 
forward way, do not limit to single out the thermodynamic properties 
of the system, but allows also to determine the dynamical evolution of 
the global quantities. It is immediate to see, for example, that, as a 
function of the virial ratio, the entropy (24) is maximized when Q = 1: 

da 2(1 - Q) 



dQ Q(2 - Q) 



(28) 



26 A very general approach to the definition of suitable microcanonical entropies 
for generic MDoF systems is discussed in (45). 

27 According to what stated above, in the derivation of the microcanonical ther- 
modynamic relations, we take into account that, for bound SGS, the energy, the 
number of bodies and the volume are not independent quantities. 
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Figure 4- Probability distribution for the values of the virial ratio, determined ac- 
cording to the definition of the microcanonical entropy, equation (24), for bound self 
gravitating systems with moderately large number of bodies: from the broader to 
the narrower distributions it is N = 10, 50, 200, 1000. It is well evident that, as the 
number of bodies increases, the overwhelming fraction of Ef is occupied by systems 
in virial equilibrium. 



That is, during the approach towards the virial equilibrium the average 
entropy of the system increases to the maximum value allowed by the 
values of the other parameters. Equivalently, it can be said that, start- 
ing from an initial state out of the virial equilibrium, entails a rapid 
relaxation to the maximum entropy state compatible with the given 
energy. That most of the volume of phase space of N-body systems is 
occupied by configurations in virial equilibrium is evident from figure 
4, where the probability distribution of the Q values is shown to be 
strongly peaked around the equilibrium value already for moderately 
large N- values. 

We think that it is very impressive to see how such a definition 
of finite time entropy is able to predict even the dynamical states 
maximizing locally (in time) the entropy. 

As a further evidence of this property, we observe that it is also 
possible to predict (and also to numerically verify), that the secular 
dynamical evolution of a gravitationally bound system leads the ratio 
ot(t) to an average steady growth, over which small quasi periodic 
modulations are imposed. The expression derived above for the mi- 
crocanonical entropy correctly describes also this secular dynamical 
evolution through locally stationary states of slowly increasing entropy! 

da N 

JT oc - > . 29 
oa a 

This means that the entropy will increase along with the growth of the 
parameter a(t) (because, supposing fixed all the other parameters, we 
have a oc Na/a). Well, it is easy to see, even analytically, that the 
development of a core-halo structure implies an increase of a(t). The 
relationship between a and a however indicates that the entropy growth 
rate, though remaining positive, slows down as a steadily increases. 

Adding a further comment to the results above, we remark that, 
while the apparent singularity of the entropy when e — > 0~ is fictitious, 
as, in this limit it is also Q — ► 2~ , the singularity for Q — > + is 
instead real, and again very enlightening on the sensible meaning of the 
microcanonical entropy here presented. Indeed the above limit describes 
a completely frozen state of the system, that is, a given configuration 
with all the bodies having zero velocity. The probability of such a state 
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is zero for any e < 0, and correctly its entropy diverges. Incidentally 
we notice that, if e — ► CP, then a frozen state become possible, and 
accordingly a can be made finite. 

We emphasize again how the present definition of microcanonical 
entropy is effective in taking into account the secular evolution of 
the system, through the full hierarchy of timescales associated to the 
singular dynamical processes occurring in the asymptotic evolution of 
SGS. It is also remarkable that, if we consider a somewhat regularized 
gravitational N-body system, where the divergence at short distance is 
removed (e.g., through a usual softening procedure, for example) and 
the unboundness of configuration space accessible is avoided, introduc- 
ing a kind of cosmic vessel, then the microcanonical entropy will still 
entail the secular evolution, but will be in this case able to single out 
an end state that is now well defined and attainable in a finite time! We 
stress that this is not at all a trivial result, because the regularization 
procedure described above does not modify the statitistical and ther- 
modynamic peculiarities of the system, which remains nonextensive 
and for which the inequivalence of ensembles remains true, as long as 
the system is gravitationally bound. In particular, such a regularized 
system will keep a negative heat capacity in a wide range of specific 
energies and the canonical ensemble description permains unable to 
grasp even the qualitative features of its Thermodynamics (44)! 



5. Comments and Conclusions. 

We presented an effective microcanonical entropy able to describe cor- 
rectly not only the thermodynamic properties os SGS, but also to out- 
line the sequence of secular evolutionary processes leading the systems 
towards states of higher and higher entropy. 

The dynamical and statistical justifications of a microcanonical de- 
scription of gravitationally bound N-body systems have been derived 
and discussed on the grounds of analytical and numerical evidences 
of their strongly chaotic properties, which allow the assumption of 
validity of the Chaotic Hypotesis of Gallavotti and Cohen. The re- 
sults obtained within the Geometrodynamical Approach to Chaos gave 
further support to the reliability of this assumption. 

The peculiar dynamical properties of gravitational interaction, the 
breakdown of the extensivity and the consequent inequivalence of clas- 
sical statistical ensembles, moreover, rule out any reliable canonical or 
grand-canonical description of the Statistical Mechanics and Thermo- 
dynamics of gravitationally bound N-body systems. 
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It can been shown that the microcanonical entropy is the sole ther- 
modynamic potential which keeps the extensivity at any finite time 
even for the non stable Newtonian interaction potential, while all ther- 
modynamic potentials defined within the canonical or grand-canonical 
ensembles lose the extensivity property and consequently lead to un- 
clear, unuseful or wrong thermodynamic relations. A simple sketch 
of the above assertions follows directly from the definitions of the 
canonical and grand-canonical potentials: 

Helmholtz free energy. 

In the canonical ensemble, given the internal energy 

U = {H)p = £ , (30) 

where average is performed on the canonical ensemble and £ is 
instead the fixed energy of the corresponding microcanonical en- 
semble 28 ; the Helmholtz free energy is related to it by the relation 

F = U-TS . (31) 

From the results obtained above, we see clearly that, for a non 
extensive system (be it a SGS or a small system) it is F oc iV 2 ; 
that is, F is not more an extensive quantity, although S it is 29 . 

Gibbs potential. 

In the grand-canonical ensemble, the appropriate potential is given 
by the Gibbs function, whose thermodynamic relationship with F 
is given by 

G = F + PV 

and, as it is easy to verify, both the terms of the sum are, for a 
bound SGS, proprotional to iV 2 . 

If a physically reasonable regularization of gravitational interaction, 
which does not modify its characteristic peculiarities, like the non 
extensivity and the non stability, is introduced, the proposed micro- 
canonical entropy keeps completely its validity and moreover becomes 
able to single out not only the hierarchy of secularly evolving meta- 
equilibrium states, in analogy with the non regolarized case, but also 

28 We are using the relations which are valid when there is the ensemble equiv- 
alence. Here we use them to further demonstrate the difficulties of canonical 
description! 

29 At this point we notice that the above derived expression for the temperature, 
equation (26), gives, T oc e. For non extensive systems this entails that also the 
product T S is proportional to iV 2 , coherently with all the above discussions. 
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to determine a non singular final equilibrium state, which was only 
asymptotically reachable prior the regolarization. 

In addition, the above regularization allows for a detailed study 
of the critical phenomena occurring in SGS (44), able to show the 
occurrence of phase transitions in gravitational N-body systems and 
also to confirm the doubtful reliability of a canonical description of the 
Statistical Mechanics and Thermodynamics of Gravity . 

The above discussion points out that the peculiarity of the collective 
properties of gravitational systems depend largely on the long range 
attractive nature of the interaction and it is almost completely indepen- 
dent on the formal singular behaviour at short distances. Taking into 
account that the latter is only a formal singularity, we can safely ignore 
it also on the grounds of the results showing its complete irrelevance for 
the statistical properties of SGS. The details on the way of coping with 
the dynamical consequences of the presence of the singularity affect 
the evolution of macroscopic quantities only on timescale which are 
longer than the binary relaxation time, in some cases even exponentially 
larger! 

5.1. Further comments and future developments. 

A further point deserves to be remarked: it is undeniable that the limits 
on which most of the standard Statistical Mechanical and Thermody- 
namic frameworks rely are of outmost relevance for the comprehension 
of the behaviour of most macroscopic systems. For the overwhelming 
largest class of MDoF systems, the thermodynamic limit allows to 
understand the basic and usual properties of macroscopic matter and 
even the rich and complex phenomenology of critical phenomena. At the 
same time, theorems assuring the existence of an invariant distribution 
are exceptionally important to single out the properties of the final 
equilibrium states of most macroscopic systems. 

Nevertheless, it is also evident from the considerations here pre- 
sented, that the operation of taking these limits (i.e., the thermody- 
namic limit and the infinite time limit), assigning to them the role 
of selection criteria in order to establish whether a Thermodynamics 
for a given class of systems can be constructed or not, can hinder 
in some cases the possibility of obtaining alternative statistical and 
thermodynamic description of MDoF systems which are certainly sin- 
gular in some respect, but possess notwithstanding well defined, though 
peculiar, thermodynamic properties, in general for many of them, or at 
least over any physically reasonable finite interval of time, as it happens 
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for SGS. 

The outcomes presented above are enlightening also for the criti- 
cal analysis of the results based in some way on the assumption of a 
canonical equilibrium of gravitational N-body systems. In particular, 
we claim that many results based on the maximization of entropy like 
quantities constructed using the single particle distribution function 
have to be reconsidered. 

In conclusion, despite the rigorous theorems, we have shown that 
a thermodynamic description of gravitational N-body systems is not 
only feasible, but also suggests a new approach to peculiar systems for 
which the standard paradigms of rigorous Statistical Mechanics do not 
apply. 

We plan to further develop the line of research here proposed in order 
to make it possible to obtain indications on the detailed evolution of 
SGS on the basis of thermodynamic arguments alone. 

In order to attain this goal it is necessary to have more informations 
on the detailed energy distribution of N-body systems and for this the 
use of extensive numerical simulations is required. 
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